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Abstract 

A family of nonequilibrium kinetic Ising models, introduced earlier, evolv- 
ing under the competing effect of spin flips at zero temperature and nearest 
neighbour random spin exchanges is further investigated here. By increasing 
the range of spin exchanges and/or their strength the nature of the phase 
transition 'Ising-to-active' becomes of (dynamic) mean-field type and a first 
order tricitical point is located at the Glauber (5 = 0) limit. Corrections to 
mean-field theory are evaluated up to sixth order in a cluster approximation 
and found to give good results concerning the phase boundary and the critical 

exponent (3 of the order parameter which is obtained as (3 ~ 1.0. 
05.70.Ln, 05.50.+q 
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I. INTRODUCTION 



Kinetic Ising models were originally intended to study relaxational processes near equilib- 
rium states . Later combinations of Glauber and Kawasaki dynamics were used success- 
fully in investigating questions about temperature driven nonequilibrium phase transitions 
In a previous paper a class of general nonequilibrium kinetic Ising models (NEKIM) 
with combined spin flip dynamics at T = and spin exchange dynamics at T = oo has been 
introduced ||, in which, for a range of parameters ( other then temperature) of the model, 
a directed-percolation-like Ising-to-active phase transition takes place. The line of phase 
transitions have been found to belong to the same universality class as the phase transitions 
occurring in the cellular automaton models introduced and investigated by Grassberger et 
al. |7|||. Numerical studies of other models showing similar type of phase transition have 
been reported recently PJTC| . 

In the present note we consider a generalized form of NEKIM by allowing for exchanges 
of arbitrary range, R. The mean-field (MF) limit of NEKIM phase transitions is reached 
when R — > oo and/or the probability of the exchanges, p ex , relative to the time scale of spin- 
flips approaches infinity. In a systematic generalized MF approach (GMF) fll|-13|, besides 
the lowest order approximation (ordinary dynamic MF, n = 1) also the second order cluster 
equations(n=2) could be solved exactly. Numerical solutions have been obtained up to sixth 
order. 

In this way we have found strong theoretical evidence for the conjecture, stemming from 
simulations ||, that the line of Ising-to-active second order phase transition points ends at 
the Glauber limit (5 C = 0, 5 being a parameter of the spin-flip transition rate of crucial 
importance here) with maximal exchange range and/or rate. It is shown here that this end 
point is of first order (tricritical point) and is described by plain MF theory. The relaxation 
time is obtained as r oc 1/| S |. GMF results show, that with increasing n the critical point 
becomes of second order and moves towards negative values of S of increasing absolute value. 



The coherent anomaly method fT^,|T5| has been used to extract the exponent (3 of the order 
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parameter from the results of GMF calculations yielding (3 = 1.0. 



II. THE MODEL 

In NEKIM we have started with the general form of the Glauber [0 spin-flip transition 
rate in one-dimension for spin Sj sitting at site i (sj = ±1): 

Wi = 172(1 + 5s^is i+1 ) ^1 - ^Si(si_i + s i+1 )\ (1) 

where 7 = tanh 2 J/kT (J denoting the coupling constant in the Ising Hamiltonian), T and 
S are further parameters which can, in general, also depend on temperature. When T = 
is taken then 7 = 1 and (1) leads to two independent rates: 

p RW = 2w m = T(l - 5), p an = w m = T(l + 5) (2) 

responsible for random walk and pairwise annihilation of kinks, respectively. T and 5 are 
constants to be varied. 

The other ingredient of NEKIM has been a spin-exchange transition rate of neighbouring 
spins ( the Kawasaki rate at T = 00): 

Wii+l = -Pex[l ~ SiS i+1 ] (3) 

where p ex is the probability of spin exchange. Prw, Pan and p ex have been chosen as nor- 
malised to unity, leading to the relation: 

p ex = 1 - 2r (4) 

The spin-exchange process induces pairwise creation of kinks in the immediate neighbour- 
hood of an existing kink: k — > 3k with probability p ex . From this process the ultimate 
development of an active phase can arise and in ref we have made the conjecture, and 
found numerical evidence for it, that prw > Pan (i-e. 5 < 0) is necessary for this to happen. 

Now we generalize the original NEKIM model by allowing the range of the spin-exchange 
to vary. Namely, eq.(3) is replaced by 
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w i,i+k — 2^ ea; [-'- _ S i S i+k\i (5) 

where % is a randomly chosen site and Sj is allowed to exchange with Si + k with probability p ex . 
Site k is again randomly chosen in the interval 1 < R , R being thus the range of exchange. 
The spin-flip part of the model will be unchanged. We have carried out numerical studies 
with this generalized model in order to locate the lines of Ising-to-active phase transitions. 
Spin-flip and spin-exchange have been applied alternatingly at each time step, the spin-flip 
part has been applied using two-sublattice updating, while making I Monte Carlo attempts 
at random ( I denotes the size of the chain) has been counted as one time-step of exchange 
updating. It is worth mentioning, that besides k —>■ 3k, also the process k — > 5k can occur 
for R > 3, and the new kink pairs are not necessarily neighbors. The character of the phase 
transition line at R > 1 is similar to that for R = 1, except that the active phase extends, 
asymptotically, down to 5 C = 0. This is illustrated in Fig.l , where besides R — 1, the case 
R = 3 is also depicted: the critical value of — 5 C is shown as a function of p ex (Fig l.a) and 
b)). Moreover, — 5 C as a function of R is also shown at constant T = .35 . The abscissa has 
been suitably chosen to squeeze the whole ( infinite ) range of R between and 1 and for 
getting phase lines of comparable size (hence the power 4 of R/(l + R) in case of Fig.l.c)). 

Besides the critical exponent a, used in identifying the phase transition points, also the 
other critical exponents characterizing the phase transition have been determined numeri- 
cally around some typical points (far from the end-points) of the phase transition lines for 
R > 1, with the same result as was obtained in || for the case R — 1: the exponents agree 
- within error - with those of Grassberger's automata ]7].|8|. 

On the phase diagrams of Fig.l two non- typical regimes can be distinguished, namely 
a). Pex ~ 0.( on Figs, l.a) and b).) Here NEKIM's behaviour is getting close to that of 
the plain spin-flip model at T = 0: the steady state is everywhere Ising-like except for the 
limit-point 5 = — 1, where p an = and the initial kink density is sustained. At this specific 
point the energy becomes conserved |1| and a phase transition takes place with a change in 



the form of the time dependence of correlations from exponential to stretched exponential. 
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This limit will not be further discussed here. 

b). Pex ~ 1-0 and/or R — ► oo (on Figs l.a)-c)). For p ex — > 1, R — 1 we have concluded 
in [6] that 5 C — > — 0, though it has been difficult to get reliable numerical estimates for the 
critical exponents of the transition due to the long transients. 

In limit b)., after each step of spin-flip ordering, maximal mixing of the neighbourhood 
of each spin follows, suggesting that a mean-field type situation takes place. It is important 
that according to eq.(4) p ex = 1 is approached together with r — > and thus p ex /F — > oo- 
(As T sets the time scale of the ordering process, its vanishing tendency enhances the effect 
of mixing). The same limit can be reached at fixed p ex by increasing R to infinity (Fig.l.c)). 

We have also checked the decrease of — 5 C with 1/R numerically at fixed T = .35,p ex = .3 
and found, over the decade of R = 4 — 40, that 

-5 C « 2.0(l/i?) 2 (6) 

reminiscent of a crossover type behaviour of equilibrium and non-equilibrium phase transi- 
tions |T7|, here with crossover exponent 1/2. It should be noted here, that to get closer to 
the expected S c mf = 0, longer chains (we used / values up to 20000) and longer runs (here 
up to t = 5 • 10 4 ) would have been necessary. The former to ensure / ^> R |]TB[ and the latter 



to overcome the long transients present at the first few decades of time steps. 

In what follows we will always refer to the MF limit in connection with p ex — > 1 (i.e. 
Pex/r — > oo), for the sake of concreteness, but keep in mind that R — > oo can play the same 
role. 

III. MEAN-FIELD THEORY AND CORRECTIONS TO MEAN-FIELD 

It is straightforward to find the MF equation for the spin-flip model alone ( at T = 0). 
By denoting the average spin density by M we get, using (1) 

dM/dt = —5TM(M 2 - 1) (7) 
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The fixed point solutions are M* = 1, —1, of which the first two are stable if 5 > 0, while 
the M* = solution is stable for 5 < 0, suggesting a (first order) order-disorder-type phase 
transition at 5 = 0. That the M* = fixed point is not an antiferromagnetic type can 
be shown by carrying out a two-sublattice MF analysis of the model [19]. First sublattice: 
odd lattice points with average magnetization M 1; second sublattice: even lattice points 
with average magnetization M 2 . The total average magnetization M = (Mi + M 2 )/2 and 
the difference of the sublattice magnetizations A = (Mi — M 2 )/2 obey the following MF 
equations: 

dM/dt = —5TM(M 2 - 1 - A 2 ) (8) 

dA/dt = ~2TA + 5TA{M 2 -1- A 2 ) (9) 

The solutions for fixed point A* ^ are: M* = 0, A* 2 = -1 - §. 

The values of 5 being restricted to 1 > 5 > — 1, the only possibility to ensure A* 2 > 
at the same time is: 5 = — 1, with A* 2 = 1. Thus we will suppose that the transition at 
S c = is of order-disorder type. A small fluctuation dM around one of the stable fixed 
points decreases as dM(xe~ l / T with r = 1/T\ 5 |. This relaxation time becomes infinite at 
the MF transition point. The corresponding critical slowing down in its vicinity explains 
the longer and longer transients observed during simulations. 

Fig. 2. serves to illustrate the MF result in comparison with results of simulation of 
NEKIM at three values of p ex . The average density of kinks at t = oo is depicted versus 5. 
MF approximation corresponds to Fig. 2. a), with a jump at 5 c mf = 0. Fig. 2.b). shows 
the behaviour of the pure spin-flip model at T = 0: the steady state is everywhere Ising-like 
(p(oo) = 0) except for 5 = — 1. Figs. 2.c).,d). and e.) are results of simulation of NEKIM 
at p ex = .02 , p ex = .9 and p ex = .98, respectively. By further decreasing (increasing) p ex , 
the NEKIM curves get closer and closer to b). ( a).), respectively. Such behaviour is in 
accordance with our expectations: it supports the MF interpretation of the high-p e:E part of 
the phase diagram. 
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We have applied the generalized mean-field calculation method, or cluster approximation 



TT| , [T2|| in the form applied for cellular automata [T3|] in order to go beyond the lowest order 
approximation shown above. 

Steady-state equations have been set up for block probabilities in n = 2... 6-th order. The 
system of GMF equations are solvable analitically also for n = 2. 

The n = 2 approximation gives the density of kinks p(oo) as : 



p(oo) 



; PRW 2 + Pan ~ PRW Pan ~ \J JqPRW + f PRW 2 Pan ~ \ PRW 3 Pan + Pan 2 ~ 2 p RW Pan 2 
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2 (| PRW 2 - PRW Pan + Pa 

(10) 

for 8 < 0. For 8 > p(oo) = 0, i.e. GMF still predicts a first order transition for 8 = 0; the 
jump in p(oo) at 8 = 0, however, decreases monotonically with decreasing T, according to 
eqs. (10) and (2). 

In order to get the n > 2 approximations, the set of GMF equations can be solved 
numerically only. We determined the solutions of the n = 3,4,5,6 approximations for the 
kink density at i)T = 0.35 (Figure 3.) and of the n = 3,4, 5 approximations at ii).r = .05 
(Figure 4.). As we can see the transition curves became continuous, with negative values 
for 5" (<5™ denotes the value of 5 in the n-th approximation for which the corresponding 
p(oo) becomes zero). Moreover, | 5™ \ increases with growing n values. As increasing n 
corresponds to decreasing mixing, i. e. decreasing p ex , the tendency shown by the above 
results is correct. 

Figs. 5. a), and 5.b). show a quantitative - though only tentative - comparison between 
the results of GMF and the simulated NEKIM phase diagrams. The obtained GMF data 
for 8™ corresponding to n = 3,4, ...,6 (r = .35) are depicted in Fig. 5 function of 

l/(n — 3), together with results of simulations. The correspondence between n and p ex has 
been chosen as the simplest conceivable one. ( Note that 8 C ^ is obtained first for n = 4.). 
The simulated phase diagram has been obtained without requiring the fulfillment of eq.(4), 
at constant T = .35. In this case the 8 C = limit, of course, is not reached and a purely 
second order phase transition line can be compared with GMF results (for n values where 



it also predicts a second order transition). Simulations for R = 3 have been found to lead 
to | 8 C | values low enough to fit GMF data. The (polynomial) extrapolation of GMF data 
to n — > oo (corresponding to p ex = , i.e. plain spin-flip), shown also in Fig. 5., could have 
been expected to approach 5 = — 1. That this is not case can be due to the circumstances 
that upon increasing n i) GMF starts here from a first-order MF phase transition, ii). which 
becomes second order and iii).GMF should end up at a pathological point, discussed shortly 
in section II., with quite unusual (and not yet cleared up) properties. 

Fig 5b). shows the n = 3,4,5 results for T = .05, which are compared now with the 
R = 1 simulation data. 

The critical exponent j3 of the order parameter has been determined processing the results 
of GMF approximation by the Coherent- Anomaly Method (CAM) [14,15]. According to 
CAM the GMF solution for kink density p at a given level of approximation - in the vicinity 
of the critical point, S c - is the product of some mean-field behavior multiplied by the 
anomaly factor a(n): 

p(n) = a(n) - , (11) 

The true critical exponent, (3, can be obtained by fitting, using the knowledge that the 
divergence of the anomaly factor scales as : 

a{n) ~ {5 n J8 c - , (12) 

as the level of approximation n goest to infinity. More precisely, for the available low level 
of approximations (n < 6), correction to scaling should also be taken into account : 

a{n) = b A^ MF + c A^ /3 « f+1 + ... , (13) 

where b and c are constants and the invariant variable 

A n = (4/5 c n ) l/2 - (m) 1/2 (14) 

is used. This new variable was introduced recently [JTSJ to avoid the ambiguity on the choice 
of the independent variable (5 «-> <5 _1 ). Using this new variable accurate estimate was given 
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for the critical exponents of the 3D Ising model [|15j and for the exponent (3 of the stochastic 
Rule 18 cellular automaton [ [20|j . 

From our GMF approximation results , as shown on Fig.3, we can use the n = 4, 5, 6 
data for the CAM analysis, while the n = 3 result can be taken to represent the lowest 
order MF approximation (with 5 C MF = 0) for a continuous transition ( no jump in p for 
n = 3). For 5 C we use the results of the polynomial extrapolation. Fig. 6. shows that in 
the n = 3 approximation the exponent (3 = 1.0064 , thus (3mf ~ 1- Graphs similar to that 
on Fig. 6. have been obtained for n = 4,5, 6, as well. Consequently, as Table | shows, the 
anomaly factor does not depend on n. This means, according to eq.(12), that the exponent 
is estimated to be equal to the 'mean-field' value (3 ~ (3mf = 1- 
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TABLES 



TABLE I. CAM calculation results 

n A™ a(n) 

4 2.49043 0.01083 

5 1.81022 0.01074 

6 1.45766 0.01079 
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IV. DISCUSSION 



The mean-field limit of the line of non-thermal phase transitions occurring in a family 
of one-dimensional kinetic Ising models has been analysed here. This line consists of second 
order Ising-to-active phase transition points which belong to the universality class found 
first by Grassberger et al The first order endpoint of this line has been found to be 

described by MF theory. Systematic generalized MF theory has been applied to treat bigger 
and bigger blocks of size n exactly in order to be able to depart from this tricritical point. 
Numerically solvable results up to n = 6 have given support of results of simulations and 
especially provided a value for the critical exponent (3 of the order parameter (density of 
kinks) (3 = 1, which is in accord with Grassberger's result : (3 = .94 ± .06. The value (3 = 1 
coincides with the MF /^-exponent for directed percolation. This is not surprising in case of 
our n = 3 result which we have used as an effective MF one for a continuous transition at 
5 = 0. As Grassberger has pointed out || in the rate (or MF) approximation there is no 
difference between models leading to the Ising-to- active transition and directed percolation. 
In this argumentation, however, the MF equation is written for the kink density (and not 
for the magnetization as in eq.(7)) and has the form: dn/dt = 2fin — 2Xn 2 , where /i and A 
are the reproduction and annihilation rates, respectively. Nevertheless, a heuristic equation 
of similar type can also be constructed in the present model using ji oc (prw — Van) as 
a rate making kink- reproduction effective (a conclusion stemming from simulations). The 
corresponding Mf critical value is then 8^f F = and (3 MF = 1. It is, however, surprising that 
higher order approximations of GMF have not given practically any deviation from the MF 
value of (3. In simulations for branching annihilating random walk with four offsprings Jensen 
@ conjectures a value (3 = 13/14 on the basis of simulation results, no theoretical motivation 
appears to exist for this value, however. To decide the question what the exact value of (3 
in this universality class is appears to be a challenging task, and probably calculating GMF 
in even higher approximations than here would be worth wile. 
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FIGURES 

FIG. 1. Phase diagram of NEKIM for 5 c (R,p ex ) is depicted for a). R = 1, b). R = 3 as a 
function of p ex (with 2T = 1 — p ex ) and for c). T = .35 (p ex = -3) as a function of (i?/l + R) 4 (full 
line). The phase boundaries have been obtained by measuring p(t), the density of kinks, starting 
from a random initial distribution and locating the phase transition points by p(t) oc t~ a with 
a = .27 ± .04. Typically the number of lattice points has been / = 2000 and averaging over 500 
independent runs has been performed. 

FIG. 2. Kink density p in the steady state as a function of 5: a). result of first order MF 
calculation with a jump at 5 = 0, reflecting the first order nature of the transition, b). the same 
quantity for the pure spin-flip model: p(oo) = p(0) at S = —1, otherwise p(oo) = 0. c).,d).,e).: 
results of simulations of NEKIM, see text. 

FIG. 3. p(oo) as a function of 5 as obtained from GMF for n = 3, .., 6 in case of T = .35 

FIG. 4. Kink density p(oo) obtained by GMF for n = 3, 4, 5 in case of ir = .05 

FIG. 5. Comparison between results of simulation and GMF: GMF results for S c are plotted as 
a function of 1/ (n — 3) , while simulation results for S c are depicted as a function of p ex at constant 
r, with R = 3 in case of T = .35 (Fig. 5. a)) and with R = 1 for T = .05 (Fig. 5.b).). 

FIG. 6. The n = 3 results of the GMF for p(oo) with T = .35, are shown on a double logarithmic 
plot as a function of | (S — 5 C ) \. The straight line gives (3 = 1.0064, which is used as a fictitious 
MF value in CAM. 
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